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During the Precambrian era. Earth’s decelerating rota¬ 
tion would have passed a 21 -hour period that would have 
been resonant with the semidiurnal atmospheric thermal 
tide. Near this point, the atmospheric torque would have 
been maximized, being comparable in magnitude but op¬ 
posite in direction to the lunar torque, halting Earth’s ro¬ 
tational deceleration, maintaining a constant day length, 
as detailed by Zahnle and Walker [1987]. We develop a 
computational model to determine necessary conditions for 
formation and breakage of this resonant effect. Our sim¬ 
ulations show the resonance to be resilient to atmospheric 
thermal noise but suggest a sudden atmospheric tempera¬ 
ture increase like the deglaciation period following a possi¬ 
ble ’’snowball Earth” near the end of the Precambrian would 
break this resonance; the Marinoan and Sturtian glaciations 
seem the most likely candidates for this event. Our model 
provides a simulated day length over time that resembles 
existing paleorotational data, though further data is needed 
to verify this hypothesis. 


an extended period of time. In our model of atmospheric 
resonance, there are effectively three outcomes for this res¬ 
onant phenomenon. 

First, the Earth could have entered a stable resonant state 
which lasted for some extended period of time before being 
interrupted by a global temperature increase, such as the 
deglaciation period following a possible ’’snowball Earth” 
event. Specifically, the Sturtian or Marinoan glaciations 
make good candidates for this breakage event. [Pierrehum- 
bert, et al. , 2011; Rooney, et al. , 2014] 

Second, the resonant stabilization could have never oc¬ 
curred, as the Q-factor of the atmosphere could have been 
too low for the magnitude of the atmospheric torque to ex¬ 
ceed that of the lunar torque, a necessary condition for a 
constant day length. 

Third, the resonance could have been of no interest, as 
atmospheric and temperature fluctuations could have been 
too high to allow a stable resonance to form for an extended 
period of time. 

We discuss the plausibility of each of these scenarios in 
greater detail below and ultimately conclude that the first 
scenario is the most likely to have occurred. 


1. Introduction 


2. Analysis of atmospheric resonance 


At some point during the Precambrian, the Earth would 
have decelerated to the point where it had a rotational pe¬ 
riod of 21 hours, which would have been resonant with the 
semidiurnal atmospheric tide, with its fundamental period 
of 10.5 hours. At this point, the atmospheric tidal torque 
would have been comparable in magnitude but opposite in 
sign to the lunar oceanic torque, which could create a sta¬ 
bilizing effect on the day length, preserving the 21 hour day 
length until the resonance was broken, as first discussed in 
Zahnle and Walker [1987]. 

The question then arises as to how the Earth broke out 
of its resonance-stabilized day length of 21 hr to progress to 
its current day length of 24hr. In general, any sufficiently 
large sudden increase in temperature will shift the resonant 
period of the atmosphere by thermal expansion (resulting in 
a change of atmospheric column height) to a shorter period, 
as described in Eigure 1, and could potentially break reso¬ 
nance, allowing for Earth to decelerate to longer day lengths. 
(Alternately, resonance could also be broken by increasing 
the lunar torque to surpass the peak atmospheric torque by 
the gradual change of the oceanic Q-factor, defined for an 
arbitrary system as 27 r • - total energy -^ though the 

very low necessary atmospheric Q factor for resonance to 
form given the current oceanic torque makes this seem a 
less likely explanation and is not explored here.) 

This study develops a model of resonance formation and 
breakage that approximately outlines the necessary condi¬ 
tions for this constant day length phenomenon to occur for 
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The details of the atmospheric tide are quite complex, 
but the essential features can be appreciated with the fol¬ 
lowing toy model of the torque. (For interested readers, a 
more complete treatment of the atmospheric tidal problem 
is given in Lindzen and Chapman [1969], most specifically 
in Section 3.5.C.) 

Given a fluid with column density po and equivalent col¬ 
umn height /lo under gravitational acceleration p, with Lamb 
waves of amplitude h ho and wavelength X ^ ho, wave 
speed of \/gho, Cartesian spatial coordinates of x, a forced 
heating term hf , and a damping factor f = (with tg de¬ 
fined as the total energy over power loss of the system, such 
that Q = oJotg), we start with the forced wave equation 
without drag (we will add this in later): 


d^h d^h d^hf 


( 1 ) 


We are interested in a heating term of the form F = 
Fq cos(2wf -I- 2kx), with Fo as the average heating per unit 
area, oj as the angular frequency, and k = at the equa¬ 

tor, with i ?0 the Earth’s equatorial radius. T^ius, for Cp as 
the specific heat at constant pressure and To as mean surface 
temperature, we have that poCpTo-^ = Fo cos( 2 cut -|- 2kx), 


Fo sin( 2 cjf -|- 2kx) 
2poCpToui 


( 2 ) 


Expressing h = Asin(2cJt -|- 2kx) and defining the equiva¬ 
lent height (which is currently 7.852km [Zahnle and Walker 
, 1987], with resonant effects occurring at about 10km) as 

ho = ) where /3 is the relevant eigenvalue to Laplace’s 

tidal equation. Using present values for ho, w, g, and k, 
we obtain (3 « 0.089, which agrees swith Lamb [1932], p. 


1 



X - 2 
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560. Near resonance, gho = 
Equation 1 that: 

A = -- 


/3fe^ 


0k 


f, so we obtain via 


loFq 


(3) 


2 poC'pTo(a;g - 

At the present day, cu < cjq, making A negative, so the 
positive peak of Asm{2ujt + 2kx) is at 2u)t + 2kx = — f. 
At noon (t = 0), this occurs spatially ata: = —^ 
or —45°. This result determines the sign of the torque, as 
the mass excess closer to the sun exists such that it is being 
pulled in the prograde rotational direction. Note that for 
period of time where the length of day is less that the reso¬ 
nant period of 21 hr, that is, for uj > wq, the resultant torque 
of A will exert a decelerating effect on the earth. However, 
at the point of resonance in question, where the lunar torque 
is cancelled by the atmospheric torque, cu < uq by a small 
factor, corresponding to a day length slightly above 21 hr. 

Addressing drag in our model, if we assume that any ex¬ 
cess velocity formed from the tidal acceleration in the atmo¬ 
sphere is quickly dissipated into the Earth through surface 
interactions with a damping factor E, and that this surface 
motion is relatively quickly dissipated into the rotational 
motion of the entire Earth, as given by Hide, et al. [1996], 
writing the dissipative Lamb wave forces, we have: 


dv dh 

iH ^ 

from which we obtain: 


dv dh dhf 


A = 


Fo {2 lv — — Lvl)-\-2iLvV) 

PoCpTo 16(a;2 _ ^2)2 + 4^2r2 • 


(5) 


In this model, the imaginary component O'(A) represents 
amplitude which would create a force with angle of ^ with 
respect to the sun, and thus does not exert any torque on 
the Earth. We need only concern ourselves with the real 
part 5R(A), then. Thus, we have: 


K(A) = 


( 6 ) 


Fo — ^o) T 

2poCpTo ' 4(a;2 -cu2)2 -^a;2r2' 

Since we know the atmospheric torque to be directly pro¬ 
portional to the atmospheric displacement A, we can use 
the fact that the present accelerative atmospheric torque, 
2.5 X 10^®Nm, is approximately that of the present de- 
celerative lunar torque, 4 x 102®Nm, as given in Lambeck 
[1980], to scale the atmospheric torque along the curve fol¬ 
lowing 5R(A), giving the total atmospheric torque Tatm{<v) as 
a function of the Earth’s rotational frequency, as detailed in 
Eigure 1. 

Eor a sufficiently high atmospheric Q, we can see that, 
starting from a short day length, as the Earth decelerates 
over time, increasing the length of day, the atmospheric 
torque increases until it eventually matches the lunar torque, 
keeping the length of day constant at this stable equilibrium. 
While there are two day lengths at which the torques are 
balanced, only the one at lower day length is stable. That 
is, infinitesimally perturbing the system about the unsta¬ 
ble equilibrium will cause the system to move away from 
equilibrium (to a longer day length). 

It should be noted that, to the authors’ knowledge, there 
is little consensus on a value of Q for the atmosphere, though 
Lindzen and Blake [1972] put Q for a period of 10.5 hours 
at about 30. Regardless, one can reasonably assume it is 
within the range of 10 — 500, so we solve the problems in 
this paper using all possible values of Q within this range. 
Ultimately, we establish a critical (relatively low) threshold, 
dependent on the lunar torque, that Q must exceed for res¬ 
onance to form - all systems with Q sufficiently past this 
threshold result in similar results. 


3. Estimation of resonance-breaking conditions 


Before solving the deglaciation timescale problem with a 
more complete computational model based on the previous 
section, we detail a less sophisticated analytical solution to 
approximate the warming timescale necessary to break res¬ 
onance. We then verify this with our computational model, 
noting that the key features of the solution are present in 
both models, albeit at different values. 

Given some increase in global temperature AT from an 
initial ’’average” temperature To, we would expect a cor¬ 
responding increase in atmospheric volume. Since the at¬ 
mosphere is horizontally constrained, this should result in 
a nearly linear increase in the column height of the atmo¬ 
sphere. This, in turn, would change the propagation speed 
of an atmospheric Kelvin wave, given by u = y/gho, and thus 
the resonance frequency of the atmosphere. A decrease in 
global temperature increases the atmospheric resonant fre¬ 
quency (thus decreasing the equilibrium length of day, shift¬ 
ing the curves to the left on Figure 1), while an increase in 
global temperature decreases this frequency. 

Suppose Earth had progressed to the stable equilibrium 
point in Figure 1. A large, fast (but non-instantaneous) 
increase in global temperature could shift this stable equi¬ 
librium point to sufficiently lower day lengths such that the 
unstable equilibrium point would be shifted past Earth’s day 
length, allowing the Earth to freely decelerate to longer day 
lengths. This change in temperature would need to be sud¬ 
den enough that the Earth’s rotation could not ’’track” this 
change. Additionally, decreasing temperature to its previ¬ 
ous value before the Earth has a chance to migrate away 
from the resonant zone could result in a recapture event. 

Let’s examine how fast this increase in temperature would 
need to occur. Eor a change in resonant frequency to pre¬ 
serve resonance throughout the duration of the change, the 
rotational frequency of the Earth must track the change in 
resonance frequency of the atmosphere, so ^ ~ 
make the simplifying approximation that the torque curve in 
Eigure 1 has zero width (removing any ’’buffer zone” about 
the stable equilibrium), then we require ^ Since the 


resonance frequency of the atmosphere is cjq = 
and ho (X T for temperature T, we can express ujo{T] 


€22 


as 


ujo = ' 2 fl 0 ° ’ "'ifh To the initial temperature. For any 

realistic changes in atmospheric temperature, ~ 1. De¬ 
noting the time over which the temperature changes by an 
amount AT as tw, we obtain: 


ghn T(t) , -;- 

duj^ _ dT{i) V _ AT gho I€22 , . 

dt dt 2 R^T(t^ AtwToR(^ 

Following the amplitude-scaling technique mentioned in 
the previous section, we know the maximum angular accel¬ 
eration of the Earth, or the fastest the Earth can ’’track” 
changes in cjq, to be: 


duj _ Tatm - Tmoon _ ( 16 -A( 

dt I® 7® 

where oJmax is the rotational frequency associated with 
the global maximum of Tatm, and jflJy is the current ro¬ 
tational frequency of the earth. Abbreviating A(ti;niax) as 
Amax and A{^^) as A 24 , we obtain that: 
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Torque (Nm) 



Day Length (hr) 


- Q = 50 

- Q = 100 

- Q = 200 

- Q = 200, AT < 0 

- Q = 200, AT > 0 

- Current Lunar Torque * -1 

- Precambrian Lunar Torque * -1 


Figure 1. Torque values for atmospheric torques assuming various Q-factors compared to lunar torque. Torques are 
scaled at the 24hr. endpoint such that they have a value that of the lunar torque, while the contour of the curve is 
determined by the A term derived in section 2. Note that the minimum value of Q required to form a resonance (the value 
such that its magnitude exceeds the lunar torque) varies linearly with the lunar torque. During the Precambrian, when 
the lunar torque was thought to be approximately a fourth of its current value [Zahnle and Walker , 1987], very low values 
of Q could have permitted resonance to form. 


ATgho/e22 _ '7”moon f -^max 

70 \ 16 • A 24 

Since Amax scales very nearly linearly with Q, we need 
only attain one value of Amax and scale it accordingly with 
Q. For example, at Q = 100, Amax ~ 27.01 • A24, and at 
Q = 200, Amax ~ 53.78 • A24. Using the value at Q = 100, 
our expression for the minimum tu, over which the tempera¬ 
ture can change by AT without breaking resonance becomes: 

( 10 ) 

2T’o'rmoon( jQQ 1) 

As shown in Figure 2, this expression indicates asymp¬ 
totes for stability-preserving {Q,tw) pairs as Q —>■ 60 and as 




Figure 2. The shortest possible stability-preserving 
warming time as a function of Q, as derived in Section 
3, with an arbitrary choice of AT = lOK. The asymptote 
arising at Q « 60 is a result of no resonance-stabilizing 
effect occurring, as the maximum value of Tatm fails to 
surpass the lunar torque. As this simple model serves 
only as an upper bound for the conditions required to 
break resonance, the asymptotic value of Q « 60 and the 
timescale will be lower in the computational model, 
but the general behavior described should still be present. 


tin —> 0. For a plausible atmospheric Q-factor of 100 and a 
temperature change of AT = lOTf, any signihcant change in 
temperature faster than on the order of 10® years will break 
resonance. 

Note that this model simply provides an upper bound 
on the time over which the temperature can be changed 
while preserving resonance and an upper bound on the min¬ 
imum threshold for Q for resonance to form. This is due to 
the approximation that the curves in Figure 1 have a half¬ 
maximum width of zero and that resonance will be broken 
Earth’s rotation at all deviates from the shifting equilib¬ 
rium. In reality, the width of the curves provides a buffer 
for the Earth; for example, very small changes of resonance 
frequency will not break resonance, as the displacement is 
not sufficient to put the Earth outside of the stable zone, 
even if the change occurs instantaneously. Thus, tempera¬ 
ture changes will need to occur over a shorter than the 
values shown in Figure 2 to actually break resonance, and 
resonance may form with lower values of Q. Additionally, 
the current ratio of atmospheric to oceanic lunar torques is 
hard-wired in this estimate, overstating the oceanic torque 
and understating the strength of the resonance. These prob¬ 
lems are more precisely addressed with our computational 
model outlined in the next section. However, this simple 
model gives us a good idea of the behavior we should ex¬ 
pect: an asymptote at low Q, t-u, scaling approximately with 
i.Q_i , and, for sufhciently large AT such that the buffer 
zone is small in comparison to the induced change in reso¬ 
nant frequency, scaling linearly with AT. 

4. Computational model 

To determine the required conditions for a temperature 
change to break resonance, a computational model was 
implemented in Python to numerically compute solutions 
to the model developed in Section 2 over the (AT,Q,t-w) 
parameter space. This generated a stability regime dia¬ 
gram depicting choices of (AT, Q, tu,) resulting in stable and 
resonance-breaking (unstable) states, shown in Figure 3. 

At the program’s core is a simulation function which iter¬ 
ates the Earth’s rotational frequency as a response to lunar 
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Minimum Warming Time: tw (Myr) 


Figure 3. The stability-instability boundary calculated along varying AT, Q, and tw values. At a fixed Q-factor and 
warming time (Qo,tu)o), the resulting critical temperature change ATcrit represents the boundary between stable and un¬ 
stable resonances: any change ATcrit + e over the same period of time two will break resonance, while any change ATcrit — e 
over two will preserve resonance. This can be visualized as a bent surface with the AT scale directed out of the page; 
stable resonances lie below this surface, while unstable resonances lie above the surface. 

Higher Q-factors permit larger temperature changes per unit time (^), as the system is more responsive to external 
torques than scenarios with lower Q. Conversely, for a fixed critical temperature change ATcrit over a period of time 
tw, smaller values of tw require the system to be more responsive to external torques to preserve resonance, requiring a 
larger Q. It should be noted that, regardless of Q and tw, there exists a nonzero minimum value of ATq required to break 
resonance (about 5.7 K in the simulation). 


and atmospheric torque as global temperature rises from 
To — AT to To (with To being an average global tempera¬ 
ture of 287K, though this precise value is unimportant) over 
an interval of tw years, simulating the warmup following a 
period of low global temperatures. 

In the absence of a reliable history of the lunar (oceanic) 
torque, the torque was simulated by taking a ’’base” value 
of lunar torque, tq, corresponding to the present-day lunar 
torque (with tq also acting as a scaling factor for Eq. 6 to 
convert atmospheric oscillation magnitude to torque) and 
scaling it by Tmoon oc tq ■ where to is the age of the 

Earth, t is the progression of the simulation (from t = 0 to 
t = to), and ro and r are the lunar orbital radii at times to 
and t, respectively. The (rather arbitrary) scaling by ^ is to 


loosely approximate the suspected time evolution of oceanic 
Qoc (and thus Tmoon) Over the history of the Earth, as Qoc is 
thought to have increased over time since the Precambrian. 
(However, as can be concluded from Figure 5, above some 
critical lower bound, the actual scaling of the lunar torque 
has little impact on the evolution of the entire system.) 

A very small step size (50yr) was used to minimize nu¬ 
merical error, particularly while simulating at very high Q 
values. The simulation function returned whether the result 
was stable (still trapped in a resonance-stabilizing region) af¬ 
ter a warmup period and a subsequent rest period to allow 
for w to settle had passed. To increase computational effi¬ 
ciency, only the stability-instability boundary (the surface 
shown in Figure 3) was solved for using a multiprocessed 
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binary search, such that the entire simulation ran in a more 
feasible 0{n^ logn) time over the parameter space. 

5. Results - tu, timescale 

A regime analysis was performed using our computational 
model to determine which combinations of atmospheric Q, 
total temperature change AT, and warming time re¬ 
sulted in a break of resonance and which preserved reso¬ 
nance, shown in Figure 3. 

As expected, for arbitrarily small t™, temperature 
changes greater than a critical threshold AT « STL will al¬ 
ways break resonance. Above AT « 20, seems to scale 
linearly with AT. The required to preserve resonance 
varies inversely with Q: with lower Q, temperature changes 
must take place over a larger period of time, as the Earth 
does not track changes in equilibrium as quickly. Addition¬ 
ally, the simulation reveals an asymptote for tw near Q = 10, 
with Q-factors below this value prohibiting resonance from 
forming in the first place. All of these behaviors are consis¬ 
tent with the model developed in Section 3. 

The overall timescale for the required to break reso¬ 
nance was smaller (by about an order of magnitude) than the 
rough estimation from Section 3: for a AT of lOK and a Q of 
100 , temperature changes occurring on a timescale shorter 
than 10^ years would be sufficient to break resonance, while 
for Q = 30, as suggested by Lindzen and Blake [1972], 
changes shorter than 30Myr break resonance, as shown in 
Figure 3. (It should be noted that while these calculated 
values of are probably correct to within half an order 
of magnitude, they are calculated with numerous assump¬ 
tions and approximations; the relative behavior is more im¬ 
portant.) Note that the break in resonance is, of course, 
conditional on the temperature staying near this increased 
temperature long enough for the Earth’s rotational velocity 
to decelerate sufficiently away from the area near resonance 
- a process which would also take on the order of 10^ years. 

Thus, our simulations indicate that, ignoring the possi¬ 
bility of recapture, had the rotational velocity and temper¬ 
ature of the Earth previously reached an equilibrium dur¬ 
ing a ’’snowball Earth”, virtually any realistic subsequent 
deglaciation period would break resonance, as discussed fur¬ 
ther in Section 8. 


I 

a. 


Figure 4. Regime analysis of sinusoidally driven atmo¬ 
spheric temperature fluctuations across (half-wave) am¬ 
plitude and frequency for an initial phase of zero. Grey 
regions indicate resonance-preserving scenarios, while 
white regions break resonance. The ’’noise” in the dia¬ 
gram, such as the small island of stability in the white re¬ 
gion is due to the fact that breakage also depends weakly 
on initial phase of the sinusoidal driver. However, phase 
was found not to change the overall shape of the curve, 
aside from small changes near the edge, so the resilience 
of the atmospheric resonance to realistic thermal noise is 
independent of phase. 


6. Results - effects of thermal noise on 
resonant stability 

In addition to a systematic global climate change follow¬ 
ing a cool-constant-warm pattern like a snowball event, the 
computational model outlined in Section 4 was also further 
developed to test the resilience of the resonance to ran¬ 
dom atmospheric thermal noise: higher-frequency fluctua¬ 
tions occurring at a variety of amplitudes. The temperature 
was driven sinusoidally across a very large range of frequen¬ 
cies and amplitudes encompassing all reasonable values for 
small-scale global temperature fluctuations. These results 
are detailed in Figure 4. It was found that, for a sinu¬ 
soidally driven global atmospheric temperature, the opti¬ 
mal fluctuation period to break resonance - that is, the fre¬ 
quency whereby the required amplitude to break resonance 
is minimized - was on the order of 10000 years. However, 
the required thermal amplitude for this value was approx¬ 
imately 20K (half-wave amplitude, so a total temperature 
oscillation of 40K), which is unrealistic, so the possibility 
of resonant break due to random thermal fluctuations was 
discarded. Further evidence for discarding this possibility 
is also provided by the results from the final figure in this 
paper. 

7. Results - simulated length of day over 
time 

Finally, we used the model from the above two sections 
to create a simulation of Earth’s length of day over its his¬ 
tory, shown in Figure 5. Given the plausibility of a snowball 
event breaking resonance, we simulated a sequence of four 
snowball events, corresponding in time and duration to an 
early Paleoproterozoic glaciation discussed in Kirschvink, et 
al. [2000] and three possible ’’snowball Earth” events dur¬ 
ing the late Precambrian: the Kaigas, the Marinoan, and 
the Sturtian glaciations. It should be mentioned that the 
time and duration (and recently, even the validity [Rooney, 
et al. , 2015]) of the Kaigas glaciation are debated; it is 
included at 920Ma and given a short duration to simu¬ 
late resonant recapture for demonstrative purposes. A va¬ 
riety of base torque values (ro)were chosen for simulation. 
Throughout the simulation, random atmospheric noise was 
also simulated as the sum of several sinusoidal drivers with 
a maximum amplitude of approximately 5K. (This temper¬ 
ature curve is arbitrary and should not be interpreted as an 
actual thermal history of the Earth’s atmosphere, which is 
not fully agreed upon; it was generated to demonstrate the 
main points of this study.) 

Existing stromatolite data as compiled in Williams [2000] 
put the point of resonant breakage near 600Ma., while the 
data points at approximately 2Ga., though even less reliable, 
could very tentatively establish a lower bound on the for¬ 
mation of this resonance. After 600Ma, stromatolite, coral, 
and bivalve data indicate that the day length increases to 
its current 24 hours day length quickly after a period of rel¬ 
atively constant day length (though paleorotational data is 
nearly absent during most of this range, only available near 
the endpoints). However, this data, particularly the early 
stromatolite data [Panella , 1972], should not be taken too 
seriously. [Zahnle and Walker , 1987] Paleontologists Scrut- 
ton [1978], and Hofmann [1973] also found this data to be 
unreliable and unsuitable for precise quantitative analysis. 
Regrettably, no significant additional data has emerged in 
the past several decades; further and more reliable data will 
be needed to test both Zahnle’s and Walker’s hypothesis and 
our developments on mechanisms of breaking resonance. 
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Figure 5. Day lengths with varying choices of tq (shown in blues and greens) and temperature values (red) over the 
lifetime of the Earth. Note that atmospheric thermal noise does not influence the day length except very near resonance, 
and that the resonant effect remains unbroken through this noise until two successive simulated snowball events at the 
end of the Precambrian 720Ma. and 640Ma., corresponding to recent estimates of the Sturtian and Marinoan glaciations. 
[Rooney, et al. , 2014] Recapture events can be seen at 870Ma following the ’’Kaigas glaciation” and, for some values of tq, 
following a Paleoproterozoic glaciation detailed in Kirschvink, et al. [2000]. Approximate empirical day length data from 
a compilation in Williams [2000] are overlayed in black (error bars included where present), and resemble the simulated 
results, though the reader should not take this data to be too reliable, particularly the data points prior to 600Ma. 


8. Conclusions 

Our model supports the first scenario presented in the 
introduction: the Earth entered a resonant state, perhaps 
at about 2Ga. before present (though this value is highly 
uncertain, as it depends on an unknown evolution of lunar 
oceanic Qoc for that epoch which is only crudely simulated 
in our model). The Earth then escaped resonance at about 
600Ma. (this value also depends on Qoc), when resonance 
was broken by a global temperature increase that is well 
explained as the deglaciation following a snowball event. 

As shown in the Figures 2 and 3, an asymptote dependent 
on lunar torque exists such that there is a critical value of 
atmospheric Q below which resonance will not form. Near 
this value, the resonance is quite unstable. Computationally, 
this asymptotic value was found to be very low: Q ~ 10 for 
the present lunar torque, and possibly lower for some smaller 
estimated Precambrian lunar torques [Zahnle and Walker , 
1987], making resonance formation likely given an estimate 
of Q « 30 from Lindzen and Blake [1972]. 

The minimum warming time required to break a reso¬ 
nance state was found to be within values that would be 
broken by a deglaciation following a snowball event; for 
most values of Q, the deglaciation period would need to 
be less than lO^yr., easily within the estimates for snow¬ 
ball deglaciations presented by Hofmann and Schrag [2002]. 
Snowball events with depressed, relatively stable tempera¬ 
tures lasting for a period of around 10^ years (also similar 
timespans as in Hofmann and Schrag [2002]) were found to 


provide sufficient time for an equilibrium of uj and cjo to be 
reached such that the subsequent deglaciation breaks res¬ 
onance, though this value also depends inversely with the 
lunar torque, which is not accurately simulated over time in 
our model. 

The mid-Precambrian was lacking in global or near-global 
glaciations, with the possible exception of the Huronian 
glaciation ca. 2.4-2.2 Ga.[Melezhik , 2006], which likely oc¬ 
curred before resonance had formed. Similar early Paleopro¬ 
terozoic glaciation was argued to be a "snowball event” in 
Kirschvink, et al. [2000]; however, the parameters described 
in that study were simulated and resulted in resonant recap¬ 
ture. The fact that there is little evidence of any potentially 
resonance-breaking glaciation for almost a billion years prior 
to the Sturtian glaciation [Rooney, et al. , 2014] lends cre¬ 
dence to the idea that the deglaciation of a ” snowball Earth” 
was the likely trigger that broke resonance after allowing it 
to persist for a length of time on the order of a billion years. 

It should be noted that while a reasonable choice of at¬ 
mospheric and lunar variables makes the scenario described 
in this study possible and likely, the paleorotational data 
available is not sufficient to confirm the hypotheses of reso¬ 
nance formation or breakage. Further data is required; it is 
our hope that this work will encourage developments in this 
area. 

Appendix A: Source code 

All of the code used in this paper is available upon request 
from the first author. 
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